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1. Introduction 

During the last dècades applications based on the most radical predictions and features 
of quantum mechanics, namely the existence of quantum superpositions and the 
entanglement, have been conceived pQ. Both of these features lie at the root of a 
new class of quantum technologies exploiting the coherent manipulation of quantum 
systems. Such process involves very delicate procedures since the inevitable interaction 
of the quantum systems with their surrounding leads to loss of information and 
causes two types of phenomena countering the coherent control of quantum systems: 
decoherence and dissipation. Decoherence, which is the loss of phase coherence between 
superpositions of quantum states, and dissipation, which is the leakage of population 
from the system to the environment, are the major enemies of quantum technologies. 
In order for quantum technologies to emerge, decoherence and dissipation need to be 
defeated or controlled, and for that they must be completely understood. 

The theory of open quantum systems deals with the microscopic description of 
quantum systems interacting with their surroundings [2]. It is of crucial importance 
to incorporate methods and tools of the theory of open quantum systems to the 
investigation of quantum technologies. Understanding and being able to describe 
the effects of the environment on quantized systems is the pre-requisite for the 
implementation of realistic applications. 

The exact analytic description of the dynamics of an open system is not an easy 
task. Commonly a series of approximations are necessary in order to obtain treatable 
equations of motions for the density matrix of the open system. Typical approximations 
are the weak coupling and the Markovian approximations. The first one assumes that 
the coupling between the system and the environment is weak while the second one 
assumes that the correlations between system and environment are negligible. 

The miniaturization process which has lead to the possibility of manipulating 
coherently quantum systems, exploiting the crucial features of quantum theory, has 
been accompanied by the need of speeding up the quantum engineering methods. 
This is because the quantum properties of the systems, such as entanglement, gets 
eventually destroyed by the interaction with the environment. The description of the 
short time dynamics of an open quantum system requires exact approaches which do 
not rely on the Markovian approximation. Recently, a non-Markovian description of 
quantum computing showing the limits of the Markovian approach has been presented 
jsJj. Moreover, for many sòlid state systems, the properties of the reservoirs are such 
that the Markovian assumption does not hold at any time. This is e.g. the case of 
atom làsers and photonic band gap materials |3j. Therefore the study of the dynamics 
of paradigmàtic open quantum systems using non-Markovian analytical approaches is 
crucial for the development of prototypes of quantum devices which are the building 
blocks of tomorrow's technologies. 

In this paper I focus on an ubiquitous open quantum system, namely the damped 
harmònic oscillator or quantum Brownian particle in a harmònic potential 0. Recently 
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an analytic approach for solving the non-Markovian master equation for this system 
via the simmetrically ordered characteristic function has been presented jü], and its 
dynamics has been studied 0. One of the new dynamical features discovered in Ref. 
[7j is the existence of a regime characterized by virtual exchanges of energy between the 
system oscillator and its surrounding. Here we focus on this regime and we present a 
method to calculate exactly the non-Markovian time evolution of the Wigner function 
of an initially squeezed state of the oscillator. It has been demonstrated in Ref. jH] 
that the regime considered in this paper may be experiment ally simulated with single 
trapped ions. It is worth recalling that, for these systems, measurements of the Wigner 
function have been performed in the experiments [Oj. 

The paper is structured as follows. In SecOwe review the properties of the system, 
introducing the master equation and its solution. Section El contains the main result of 
the paper, namely the derivation of the exact dynamics of the Wigner function for an 
initially squeezed state. Finally Sec. |U presents conclusions. 



2. The system and the master equation 



The reduced density matrix p$ describing the dynamics of the quantum Brownian 
particle obeys, in the secular approximation, to the following non-Markovian master 
equation [Tü] 

dp s A(í)+ 7 (í) 



dt 



A(t)- 7 (í) 



a)aps 



ao) ps 



2ap s a [ + psa)a 
2a) psa + psaa) 



(1) 



where we have introduced the bosonic annihilation and creation operators a and a) , 
respectively. The coefficients A(t) and 7 (i) are the diffusion and dissipation coefficients, 
and they are defined in terms of the noise and dissipation kernels, respectively j7]. The 
form of Eq. (0) is similar to the Lindblad form [llj, with the only difference that the 
coefficients appearing in the master equation are here time dependent. We say that this 
master equation is of Lindblad-type when the coefficients A(í) ± 7 (í) are positive at all 
times [T2]. Note, however, that Lindblad type master equations, contrarily to master 
equations of Lindblad form, in general do not satisfy the semigroup property. 

For a high temperature reservoir with Ohmic spectral density the analytic 
expression of the diffusion and dissipation coefficients is the following 

r.,2 



A(t) = 2g 2 kT 



1 + r : 
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1 + r 2 

with g system-reservoir coupling constant, r = uj c /lü , oj c cutoff frequency of the 
spectrum of the reservoir, uq frequency of the system oscillator, and k Boltzmann 
constant. 
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The solution for the density matrix of the system may be written in terms of the 
simmetrically ordered characteristic function Xt(0 & t time t, also known as quantum 
characteristic function (QCF), defined through the equation [T3j 

Ps(t) = ^ fxt(Oe^-^d% (4) 

Z7T J 

Using Eq. (J1J) one may solve the master equation given by Eq. (pQ) by means of a method 
based on the algebraic properties of the time evolution superoperator. Indeed we have 
proved that the form of the time evolution superoperator may be factorized into three 
parts allowing the derivation of an exact analytic solution jü]. The QCF of the system 
considered here is given by 



Xt(0 = e~^\o [e-^/ V*nj , (5) 

with Xo QCF of the initial state of the system. The quantities Ar(£) and T(t) appearing 
in Eq. (JSJ are defined in terms of the diffusion and dissipation coefficients A(í) and j(t) 
respectively as follows 

r(t) =2Í\(t 1 )dt u (6) 

Jo 

A r (í) = e" r W [* ^WAfàdtL (7) 



It is well known that non-Markovian features usually occur in the dynamics for 
times t <Cíj{ = 2tc/uj c , with tr reservoir correlation time. For quantum optical systems, 
since üü c 3> ojq, and üüq is an optical frequency, deviations from the Markovian dynamics 
appear for very short times, usually too short to be measured in the experiments. With 
trapped ion systems ujq ~ 10 7 Hz, and therefore non-Markovian features would appear 
for times t -C 0.1/ís. This is the reason why typical trapped ion experiments, wherein 
the time scales go from 1 to 100/ís, show a Markovian behaviour. In the trapped ion 
context, however, it is possible to engineer an artificial reservoir wherein parameters as 
the system- reservoir coupling and the cutoff frequency can be controlled jHj. With an 
appropriate engineered reservoir one may force non-Markovian features to appear by 
'detuning' the trap frequency from the reservoir spectral density [8j. This corresponds, 
for example, to the case in which r = uj c /üúq <C 1, e.g. r = 0.1. In this case the 
reservoir correlation time is bigger than the period of oscillation of the ion, = 27t/ís, 
and therefore the non-Markovian features show up in the time evolution and can be 
measured. At the moment, in fact, this would require a slight modification of the 
current experimental apparatus used in the group of David Wineland at NIST since a 
switch should be positioned between the filters, which are inside the vacuum chamber, 
and the trap elèctrode [T5| . 

In the following we will pay special attention to the experiment ally realizable regime 
r< 1. It is easy to verify that for r < 1 the coefficients A(í) ± j(t) appearing in the 
master equation, with A(í) and 7(í) given by Eqs. (0)-©; take negative vàlues [see, 
e.g. Ref. [12] ]• in this case, therefore, the master equation is of non-Lindblad-type. 
It is worth underlining that, in this regime and for the system here considered, the 
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Figure 1. Time evolution of the average quantum number (n) vs the dimensionless 
time r = u) c t for an initial state having excitation number (n)(t — 0) = 3. The 
parameters used are the following r = 0.05, g = 0.1, lo c /2itKT = 3 • 10~ 5 



density matrix is completely positive at all times t while the semigroup property of 
the generator of the dynamics is clearly violated 0. Under these conditions virtual 
exchanges of energy between system and reservoir appear due to the system-reservoir 
correlations. These virtual processes show up in the dynamics of the mean energy of 
the system in the form of oscillations (see Fig. QJ), and they are a clear signature of the 
break down of the semigroup property. The virtual exchanges of excitations characterize 
the regime r < 1 since in this regime the period T = 2tt/ujq of the system oscillator 
is smaller than the reservoir correlation time tr = 2it/uj c . Due to the time-energy 
uncertainty principle, for times t < T, a virtual process consisting in the absorption of 
a quantum of energy Huq from a high temperature reservoir and conseqüent re-emission 
of the same quantum of energy may occur. In Fig. ^we see the oscillation in the mean 
number of quanta (n) of the system oscillator for an initial state having (n)(t = 0) = 3, 
during the time interval < r < 1, with r = oü c t, for r = 0.05. Note that the period of 
the system oscillator, in the same time units, is u c T = 2irr ~ 0.3141, which coincides 
exactly with the period of the oscillations displayed in FigO 

In the following section we will see how virtual processes influence the dynamics of 
the Wigner function. 

3. The time evolution of the Wigner function 

The Wigner function is the Fourier transform of the simmetrically ordered characteristic 
function 




(8) 
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Inserting Eq.flSJ) into the previous equation one obtains 

?2 



W,(a) 



-I d^e- A ^ 

TC Z J-oo 

Inserting the inverse Fourier transform of Eq. 

1 
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into Eq.(JSJ) gives 
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A r (í) 
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with 



b{t,a,a ) = a-a e- m/2 e iul0t . 



(10) 
(11) 



In the derivation of Eq. (jlOj) we have used the property that the Fourier transform of 
a Gaussian is a Gaussian. The quantity ^(«(ao) is the propagator which, for t — > 0, 
tends to the delta function ó(a — ao). 

If the state of the input field is a coherent state |a ), then the integral appearing 
in Eq. (fTÜ|) may be easily calculated since the Wigner function of the initial coherent 
state is a Gaussian. In this case the Wigner function at time t reads as follows 

2" 



Wt(a) 



7r[A r (í) + l/2] 



cxp 



a^e 



-r(t)/2 p -iuiot 



oc 



A r (t) + 1/2 



(12) 



Having in mind this equation, and with the help of Eqs. ©-© and (j2J)-©, it 
is possible to show that the system-reservoir interaction spreads the initial Wigner 
function. Breathing of the Wigner function, that is the oscillation in its spread, appears 
in correspondence of the virtual processes. This is a dynamical feature which is absent 
both in the Markovian dynamics of the damped harmònic oscillator and in the Lindblad- 
type non-Markovian dynamics. Indeed, in both the previous regimes, the spread in 
the Wigner function simply increases, linearly in time in the Markovian case, and 
quadratically in time in the non-Markovian Lindblad-type case. 

We now consider the case of an initially squeezed state. The initial QCF for 
squeezed coherent state is 



Xo(0 = exp 



(13) 



Here C s = cosh(s) and S s = sinh(s), ao is the displacement of the initial state of the 
oscillator and z = se~ 1 ^ is the squeezing argument. For this initial condition the explicit 
integration appearing in Eq. (jlUj) is more difficult to carry out. In the case of an initial 
squeezed vacuum state («o = 0) with squeezing angle = 0, the Wigner function at 
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time t has the form 



1 r°° 
W t (a) = - 

71 J — oo 



e -A r (í)|Ç| 2 e (ar-«*Ç) 



exp 



-r(t)u -w t 



(14) 



Looking at the previous equation one realizes that what we need to calculate is the 
Fourier transform of the product of two functions 

1 



where 



W,a) 



7T 



^e/i(0/ 2 (0e 



(a£*-a*£) 



-Ap(í)|í| 5 



/ 2 (0 = exp 



(15) 

(16) 
(17) 



If we denote the Fourier transform JF of the functions and ,/ 2 (£) as follows, 

F\ = J-'(fi) and F 2 = J r (f2), we can recast W t in the form 

W t {a) =F 1 *F 2 (18) 

where * denotes the convolution. We now extend the method developed in Ref. ^üj for 
the calculation of the Wigner function of an initial squeezed state to the non-Markovian 
case we are dealing with in this paper. Carrying out the calculations, in a frame rotating 
with u , one obtains 

-2a 2 



W t (a) = Mexp 
x exp 



2A r (í) + e- r «(C 2s + S 2s )- 1 



2A r (í) + e- r W(C 2s -5 2s )- 1 



Mexp 



—at 



+ 



-a: 



(19) 



(Ax) 2 (í) {Ayf{t)_ 

Here, a x and a y are the real and imaginary parts of a, M is a time dependent 
normalization constant, and 

-r(í) -2b 

(20) 
(21) 

(22) 
(23) 



2 

e -T{t) e 2s 



(Ax) 2 (t) = A r (í) + 

(Ay) 2 (í) = A r (í) + 

are the variances of the dimensionless quadratures 

x = —= (a + 
v/2 V 



y 



—i 

71 



(a — a)^j 



At t = 0, (Ax) 2 (0) = e~ 2s /2 = a 2 /2 and (Ay) 2 (0) = e 2s /2 = 1/2ct 2 . We consider an 
initial state displaying squeezing in the x quadrature, i.e., in the trapped ion case, a 
state squeezed along the x direction of the trap. This corresponds to o < 1 since the 
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Figure 2. Time evolution of (Ax) 2 vs the dimensionless time r = w c í. The parameters 
used are the following a 2 = 0.1, r = 0.05, g = 0.1, u c /2wKT = 3 • 10~ 5 . The insets at 
the bottom of the figure show the contour plot of the Wigner function at times r = 0, 
t = 0.15, r = 0.3, r = 0.45. 

x-squeezing is characterized by (Ax) 2 < 0.5. In Fig. |2]we plot the time evolution of 
(Arc) 2 as a function of the dimensionless time r = u c t, for an initial squeezed state 
with cr 2 = 0.1, together with the contour plot of the Wigner function at four different 
time instants, r = 0, r = 0.15, r = 0.3 and r = 0.45. The graphic shows that, due to 
the interaction with a high temperature reservoir, the variance of the x dimensionless 
operator, defined by Eq. (j22j) . increases from its initial value (Ax) 2 (0) = 0.05. For short 
times, however, the dynamics is governed by oscillations between a squeezed state and 
a non-squeezed state due to the occurrence of virtual exchanges of excitations between 
the system oscillator and the reservoir. As clearly shown in the figure, at r ~ and 
r ~ 0.15 the system is in a squeezed state [(Ax) 2 (0) < 0.5] while at r ~ 0.3 and r ~ 0.45 
the system is in a non-squeezed state. By using the calculation we have carried out in 
this paper for the time evolution of an initially squeezed state, see Eq. (JTÏÏJ), we can also 
plot the contours of the Wigner function at these time instants [see insets at the bottom 
of of Fig. ©]. More in general Eq. (0, together with Eqs. ©-© and teli 
us that the Wigner function of an initially squeezed state, in the phase space, spread 
along both the a x and the a y axes and eventually approaches the Wigner function of 
a thermal state. For times shorter than the reservoir correlation time, however, in the 
non-Lindblad regime both (Ax) 2 and {Ay) 2 oscillate, so that if the system is initially 
in a squeezed state along the x (y) direction then this may result in x (y) squeezing 
oscillations. 

4. Conclusions 

In this paper I have studied the non-Markovian dynamics of a quantum harmònic 
oscillator in phase space. I have presented a method to calculate the Wigner function 
of an initially squeezed state of motion for short times, when the system-reservoir 
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correlations dominate the dynamics. The solution obtained is used to investigate 
the non-Lindblad regime, characterized by virtual exchanges of energy between the 
system and the reservoir. It is shown that the virtual processes are responsible for 
the oscillations between a squeezed and a non-squeezed state of the system. It is 
worth recalling that the non-Lindblad regime is experiment ally realizable with artificial 
reservoirs in the trapped ion context. In this context, moreover, the Wigner function 
is experimentally measurable. Therefore, by measuring the Wigner function, it is 
possible to reveal the virtual processes characterizing the non-Markovian non-Lindblad 
dynamics, and in particular the squeezing-nonsqueezing oscillations. The result here 
presented contributes to the understanding of the short time dynamics of an ubiquitous 
model of an open quantum system, namely the damped harmònic oscillator, and hence 
it is also of potential interest for quantum technologies. In fact in many solid-state 
physical systems promising for building prototypes of quantum devises the Markovian 
approximation does not hold. Moreover, if one is interested to very short system 
dynamics one needs to abandon the Markovian assumption. The result has also 
implications in the study of fundamental issues of quantum measurement theory since 
there are indications that the non-Lindblad regime may correspond to a situation in 
which no measurement scheme interpretat ion of the non-Markovian stochastic process 
exists. 
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